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MACHINE CAPABILITY VERIFICATION AND DIAGNOSTICS (CAP/DIA) 
SYSTEM, METHOD AND COMPUTER PROGRAM PRODUCT 

CROSS-REFERENCE TO RELATED APPLICATIONS 
[0001] This application is a continuation of U.S. Application Serial No. 
10/385,342, filed March 10, 2003, which claims the benefit of U.S. Provisional 
Application No. 60/444,225, filed January 31, 2003. 

FIELD OF THE INVENTION 
[0002] The present invention generally relates to systems, methods 
and computer program products for use in conjunction with computer 
numerically controlled (CNC) machines and coordinate measuring machines 
(CMM). More particularly, the present invention relates to a method for 
verifying the accuracy of a CNC machine or a CMM. Inherent in this 
methodology is the ability to 1) diagnose mechanical problems with the CNC 
machine or CMM that effect the repeatability and accuracy of such machine 
tools; 2) identify degradation in the performance of such machine tools that is 
indicative of mechanical failure; and 3) improve the accuracy and repeatability 
of such machine tools in certain situations. 

BACKGROUND OF THE INVENTION 
[0003] CNC and CMM machine tools are well known in the art and 

employed in virtually every field of manufacturing for tasks that include, for 
example, material cutting, material removal and inspection. Typically, these 
machine tools have one or more movable elements, such as a work table or 
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tool spindle, with each movable element being movable along one or more 
orthogonal axes and/or about a rotational axis. 

[0004] As is well known in the art, the quality of the output from a CNC 
machine or CMM is a highly complex function that is based in part on the 
accuracy and repeatability with which the movable elements of the machine 
tool are movably positioned. Those of even the most basic skill in the art wil| 
readily appreciate that the quality of the output is influenced by diverse other 
factors, such as the accuracy and repeatability of the fixturing that is 
employed to hold a tool bit or workpiece, but these other factors are outside 
the scope of the present application and as such, will not be discussed in 
detail herein. 

[0005] One way in which accuracy and repeatability may be obtained is 
through relatively tight tolerancing of the various critical machine tool 
elements, such as ball screws and ways, that control the movement and/or 
positioning of a machine toofs movable element. Issues of wear aside, one 
drawback with this proposition concerns its economics; generally speaking, 
tolerances may only be tightened to a degree without unreasonable cost 
penalty, but thereafter substantial costs are incurred so that progression (i.e., 
the tightening of the tolerances) beyond a point is cost prohibitive. With that 
in mind, the various critical machine tool elements of most if not all CNC 
machines and CMM's used in manufacturing are not capable of positioning a 
machine tool's movable elements in an extremely accurate manner (i.e., on 
the order of one ten-thousandth of an inch) without some sort of electronic 
compensation. 
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[0006] Typically, a CNC machine or CMM utilizes a controller that 
includes two sets of registers for each axis on or about which a movable 
element may be moved. The first set of registers usually employs a single 
register for use in compensating for reversal error, which is also known as lost 
motion or backlash. Reversal error occurs when the movable element is 
moved in a first direction along or about the axis and thereafter moved in a 
second direction opposite the first. The second set is a series of registers that 
correspond to the position of the movable element along or about the axis at 
various predetermined points. The registers of the second set are intended to 
correct for errors in the positioning of the movable element at the various 
predetermined points over its full range of motion along or about the axis. 
[0007] Various methodologies have been developed to quantify the 
• accuracy and repeatability of a CNC machine or CMM. Generally, these 
standards have been developed by professional societies, such as the 
American Society of Mechanical Engineers (ASME), and standards 
organizations, such as the International Standard Organization (ISO), the 
VDI/DGQ German Standard Organization, the Association for Manufacturing 
Technology (AMT) (formerly the National Machine Tool Builders Association 
(NMTBA)) and the National Institute of Standards and Technology (NIST) 
(formerly the American National Standard Institute (ANSI)). Software tools 
that provide information on the formats of such standards are well known and 
commercially available. 

[0008] While the details of these methodologies vary somewhat, they 
share analysis techniques wherein error is quantified as a band or range and 
wherein the values for the accuracy and repeatability of the machine tool are 
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based on the width of the error band. I have found, however, that values for 
accuracy and repeatability calculated in this manner are too conservative to 
be of any meaningful value. 

[0009] The plot illustrated in Figure 1 shows the results of an accuracy 
and repeatability analysis (hereinafter capability study) on a CNC machine 
using one commercially available software tool. In this study, the error band 
is based on a spread of six standard deviations (i.e., ± 3 standard deviations) 
about the mean of each test point, with the purported accuracy and reliability 
being based on the maximum width of the error band. In this example, the 
axis of the machine tool is stated to have an accuracy of 0.003850 inch and a 
repeatability of 0.003650 inch. Stated another way, the results of this 
capability study imply that if the movable element of the machine tool were 
told or commanded to move along the axis to a given point (X), the actual 
position of the movable element would be within a band of 0.003850 inch (i.e., 
X ± 0.001925 inch) and that successive repositioning of the movable element 
to the given point would place the movable element within a band of 0.003650 
inch (i.e., [X ± 0.001925 inch] ± 0.001825 inch). Accordingly, one would 
ordinarily conclude that this axis of the CNC machine would only be useful in 
manufacturing if the tolerances on the part to be fabricated were larger than 
about ±0.004 inch. 

[0010] As the data in the plot clearly illustrates, however, this 

conclusion is too conservative. Let us assume, for example, that we were to 
use this CNC machine to fabricate a part that is only two inches in length and 
that we position the part at the far left of the movable element. Since all of the 
data points from the capability study fall within an approximately ±0.0005 inch 
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band in this region, the performance of this axis of the CNC machine under 
these circumstances will be considerably better than ±0.004 inch. 
[0011] In view of the drawbacks of the known methodologies for 
quantifying the accuracy and repeatability of a CNC machine or CMM, as well 
as the common tendency of these methodologies to ignore the calibration 
process (i.e., the initializing or setting of the first and second registers which is 
discussed in detail, below), I had found no known method for programming or 
initializing the first and second sets of registers that produced accurate and 
repeatable results prior to developing my invention. The methodology most 
commonly utilized consisted of measuring the error in the positioning of a 
movable element along or about an axis at each of location designated in the 
second set of registers. This methodology appears to make use of the 
mechanic's conventional wisdom wherein the mechanic programs a given 
register in the second set to correct for measured error at a corresponding 
location of the movable element. For example, let us assume that the 
movable element is a rotary spindle that is movable in a generally vertical 
direction (Z axis) and that a given register in the second set corresponds to a 
positioning of the spindle at a location 0.50000 inch below a predetermined 
datum. When the spindle is told or commanded to move 0.50000 inch below 
the datum, the spindle does not move to that exact position - there is some 
error (let us assume for this example that the spindle moves to a point 
0.50100 inch below the datum). Using his conventional wisdom, the 
mechanic corrects for the error in positioning by entering a corresponding 
offset into the given one of the second registers. In the example provided, the 
offset would be -0.00100 inch. 
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[0012] As this process is repeated for each of the registers in the 
second set, and as there are often times 50 or more (often times hundreds) 
registers in the second set for a single axis, the process is rather time 
consuming and costly (both in terms of the cost of the mechanic's time as well 
as in machine down-time where it cannot be used in a productive manner). 
Considering that many CNC machines and CMM's used in manufacturing are 
of the type having several movable elements that are collectively movable 
along or about several axes, the amount of time to calibrate a single machine 
tool can be upwards of several days. 

[0013] With the machine tool thus calibrated, one would think that the 
movable elements of the machine tool would be capable of being moved with 
great accuracy and repeatability. Mechanics, however, have found that this is 
not the case; rather, they have observed that changes to the values in the 
registers of the second set, which are often times significant in magnitude, are 
nearly always necessary on each successive test to eliminate the positioning 
error at each given point. 

[0014] From the foregoing, it is apparent that there remains a need in 

the art for an improved methodology by which a CNC machine or CMM may 
be calibrated and its accuracy and repeatability more accurately determined. 
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SUMMARY OF THE INVENTION 
[0015] In one preferred form, the present invention provides a method 

for verifying the accuracy of a CNC machine or a CMM. The methodology 
segregates the measured error into an assignable cause portion and a 
common cause (or random error) portion. The methodology may be 
employed to affect the calibration of the machine tool to factor out the mean 
value of the assignable cause portion. Inherent in this methodology is the 
ability to diagnose mechanical problems with the CNC machine or CMM that 
effect the repeatability and accuracy of such machine tools; identify 
degradation in the performance of such machine tools that is indicative of 
mechanical failure; and improve the accuracy and repeatability of such 
machine tools in certain situations. 

[0016] Further areas of applicability of the present invention will 

become apparent from the detailed description provided hereinafter. It should 
be understood that the detailed description and specific examples, while 
indicating the preferred embodiment of the invention, are intended for 
purposes of illustration only and are not intended to limit the scope of the 
invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0017] Additional advantages and features of the present invention will 

become apparent from the subsequent description and the appended claims, 
taken in conjunction with the accompanying drawings, wherein: 
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[0018] Figure 1 is a plot illustrating several sets of error data for an axis 

of a machine tool as well as a conventionally calculated error band about the 
sets of error data; 

[0019] Figure 2 is a perspective view of an exemplary gauging system 
for performing the methodology of the present invention, the gauging system 
being shown in operative association with a conventional CNC machine tool; 
[0020] Figure 3 is a plot of a single set of error data which graphically 

illustrates a known method (referred to herein as "conventional wisdom") for 
electronically compensating for error along an axis of a machine tool; 
[0021] Figure 4 is a plot similar to that of Figure 1 , but which graphically 

illustrates two types of assignable cause error; 

[0022] Figure 5 is a schematic illustration in flowchart form of the 

methodology of the present invention; and 

[0023] Figure 6 is a partial schematic illustration in flowchart form of a 

second embodiment of the methodology of the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0024] With reference to Figure 2 of the drawings, an exemplary CNC 

machine is generally indicated by reference numeral 8. The CNC machine 8 
is illustrated to be a conventional vertical machining center of the type that is 
well known in the art. The CNC machine 8 includes a base 10, a first 
movable element or tool head 12, which is movable in two orthogonal 
directions relative to the base 10 as identified by axes X and Y, and a second 
movable element or spindle 14, which is moveable relative to the base 10 
along an axis Z, and a controller 16. Although the CNC machine illustrated 
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and described herein is a vertical machining center, it will be understood that 
the teachings of the present invention have applicability to all types of CNC 
machine tools, including but not limited to CMM's, horizontal mills, vertical 
mills, gantry mills, post mills, spar mills, wire electro-discharge machining 
(EDM) machines, vertical EDM's, water and plasma jets, routers, saws, etc. 
Furthermore, the particular CNC machine illustrated and described herein 
includes movable elements that are movable about one or more linear axes, 
but those skilled in the art will appreciate that the teachings of the present 
invention are also applicable to movable elements that are movable about a 
rotary axis. 

[0025] The controller 16 is conventionally operable for controlling the 
several drive mechanisms (not shown) that are associated with the tool head 
12 and spindle 14, as well as determining the locations of the tool head 12 
and spindle 14. More specifically, the controller 16 is operable for controlling 
the movement of each movable element along its respective axes, as well as 
continuously determining the location of each movable element along its 
associated axis or axes. The controller 16 includes a programmable memory 
18 having a first set of registers 20 and a second set of registers 22 for each 
of the axes X, Y and Z. The first and second sets of registers 20 and 22 are 
utilized to electronically compensate for error in the positioning of the movable 
elements along or about an associated axis. The first set of registers 20 
includes a single register 24 that is employed to compensate for reversal error 
on an associated axis. The second set of registers 22 includes a plurality of 
registers 26 x> with each of the registers 26 x being associated with a 
predetermined point (x) on the associated axis (or about the associated axis 
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in the case where the movable element is rotatable). In the particular 
example provided, the distance (d) between each register 26 is 0.4096 inch 
so that register 26i is associated with a distance of 0 inches from a datum 
along axis X, while registers 26 2 , 26 3 and 26 n are associated with distances 
0.8192 inch, 1.2288 inches and (n x 0.4096) inches, respectively, from the 
datum along axis X. 

[0026] A gauging system 30 is employed in conjunction with a movable 
element to sense the actual position (AP r j) of the movable element along or 
about an associated axis [the letter (r) in the subscript is utilized throughout 
this document to identify a particular run or predetermined distance over 
which data is collected, whereas the letter (i) is employed to identify a 
particular location in the predetermined distance]. In the particular example 
provided, the gauging system 30 is illustrated to include a target 32, a 
reflector 34, a laser interferometer 36 and a computer system 38, all of which 
are conventional and well known in the art. The target 32 is configured to be 
fixedly but removably coupled to a movable element (e.g., the tool head 12) 
and includes a reflective surface 40. The laser interferometer 38 is positioned 
in a convenient manner, such as perpendicular to the axis of the machine tool 
that is to be checked and at a location that is further along the axis that the 
movable element can travel. The reflector 34 is positioned in a manner that is 
aligned with the laser interferometer 36 and also along the axis in alignment 
with the reflective surface 40. As will be apparent to those of ordinary skill in 
the art, the laser interferometer 36 generates a laser beam 36a that is 
directed toward the reflector 34. The reflector 34 re-directs the laser beam 
36a to the target 32. Thereafter, the reflective surface 40 reflects the laser 
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beam 36a back to reflector 34 and then to the laser interferometer 36. The 
laser interferometer 36 utilizes the length of time that is necessary for the 
laser beam 36a to travel the route to and back from the target 32 to determine 
the position of the target 32 along the axis relative to a predetermined datum 
D. For purposes of simplifying the explanation herein, the datum D will be set 
to a value that is associated with the location of the movable element along 
the axis (e.g., if the controller 16 has moved the movable element to a 
location that corresponds to a distance of zero (0) inches along the axis, the 
datum D will be set to zero (0) in this example). 

[0027] In the example provided, the computer 38, such as a laser 

interferometer controller, is employed to determine and record the error (E rti ) 
between the actual position (AP r ,j) of the movable element along the axis (i.e., 
the position of the target 32 along the axis relative to the datum D as 
determined by the laser interferometer 36) and the desired position (DP r|j ) of 
the movable element along the axis (i.e., a point to which controller 16 has 
presumably moved the movable element) [i.e., E r .i = AP r j - DP r ,j], as well as to 
analyze the error according to the methodology of the present invention, 
which is set forth in detail, below. 

[0028] Returning to Figure 1 and with additional reference to Figure 3, I 

have found that the known methods for calibrating the axis of a CNC machine, 
such as the X axis of CNC machine 8, fail because they attempt to 
electronically cancel out all error without regard as to where or how the error 
is being generated. Figure 3 schematically illustrates the electronic 
compensation that is typically employed by mechanics when they calibrate a 
machine tool; each vertical arrow 50 represents an electronic compensation at 
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a given point along or about a particular axis, with the magnitude of the 
compensation being equal to the measured error in the positioning of the 
movable element. 

[0029] In my research, I have discovered that the aforementioned 

procedure does not produce results that are repeatable or accurate because it 
1) attempts to factor out (i.e., eliminate) all variation; and 2) fails to 
differentiate between common cause variation and assignable cause 
variation. With regard to the first point, I have found that attempts to factor out 
or eliminate all variation from a process merely result in driving the process 
out of statistical control. As such, my methodology does not attempt to 
eliminate all variation, but rather corrects for the mean value of a certain type 
of variation. 

[0030] Regarding the second point, those skilled in the art will 

appreciate that common cause variation describes random variability that is 
inherent in the process and affects all values (i.e., the absolute positioning of 
the movable element), whereas assignable cause variation is due to specific 
circumstances that can be accounted for, and whereas total variation includes 
both common cause variation and assignable cause variation. In my analysis 
and research, I have discovered two sources of assignable cause variation 
that had not been accounted for in the known methodologies for calibrating 
CNC machines and CMM's. The first source of assignable cause variation I 
have termed "reversal error" and the second source of assignable cause 
variation I have termed "cumulative error". In Figure 4, reversal error is 
generally indicated by reference letter E RE v, while cumulative error is indicated 
by the line identified by reference letter E c . 



HD&P Ref. 7784-000547/COA 
Boeing Ref. No. 00-220B (014666)) 

[0031] Generally speaking, an average value of E RE v (Erev) is utilized 
in the register 24 to compensate for reversal error, while the slope (m) of the 
line E c is utilized to compensate for cumulative error in each of the registers 
26 x . 

[0032] With reference to Figure 5, the methodology of the present 
invention will be described in detail. Although the methodology will be 
described in the context of a linear axis (i.e., axis X on CNC machine 8), those 
skilled in the art will appreciate that the teachings of the present invention 
have applicability to an axis about which a movable element may be rotated. 
[0033] The methodology begins at bubble 100 and progresses to block 
104 where a predetermined distance over which the movable element may be 
moved along or about the axis is selected and segregated into a plurality 
(N+1) of equally spaced apart axis increment points Al each of which 
separated by an interval of length or distance L. In the example provided, the 
distance is a length that corresponds to approximately the full length of travel 
over which the movable element traverses along the X axis. 
[0034] As the Central Limit Theorem is employed at a later point in the 
analysis, the quantity of axis increments Al should be relatively large so as to 
satisfy generally accepted statistical principles. I suggest a length (L) that 
would result in approximately 20 to 30 equally spaced-apart axis increment 
points Al. In the particular example provided, the predetermined distance is 
120 inches (i.e., the tool head 12 may be moved 120 inches along the X axis) 
and as such, I chose a length (L) of 6 inches to provide twenty-one equally 
spaced axis increment points Al (i.e, axis increments Al 0 through AI 2 o). 
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[0035] The methodology next progresses to block 106 where the 
values of (r) and (i) are set to one (1) and zero (0), respectively. Thereafter, 
the methodology progresses to block 108 wherein the controller 16 of the 
CNC machine 8 is employed to cause the movable element (i.e., the tool head 
12) to move in a first direction to a desired position DP rj that corresponds to 
an associated one of the axis increment points Alj. The methodology then 
proceeds to block 112 where the positioning error E r j is determined and 
recorded. As noted above, the positioning error E r j is the deviation between 
the desired position DP r>i of the movable element and the actual position AP rJ 
of the movable element. Thereafter, the methodology proceeds to decision 
block 116. 

[0036] In decision block 116 the methodology determines whether the 

value of (i) is equal to the value of (N) to thereby determine whether error data 
has been collected for an entire run (r). If the value of (i) is not equal to the 
value of (N), then the methodology proceeds to block 118 where the value of 
(i) is incremented by one (i.e., i = i + 1). Thereafter, the methodology loops 
back to block 108 so that the controller 16 may increment the tool head 12 to 
a location associated with the next axis increment Al i+ i. 

[0037] Returning to decision block 1 16, if the value of (i) is equal to the 

value of (N), then the methodology proceeds to block 120 wherein the value 
of (r) is incremented by one (i.e., r = r + 1). The methodology then proceeds 
to block 122. 

[0038] In block 122, the methodology increments the run (r) (i.e., r = 

r+1) and the controller 16 of the CNC machine 8 is set or switched to cause 
the movable element (i.e., the tool head 12) to move in a second direction 
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opposite the first direction to a desired location that corresponds to an 
associated one of the axis increments Alj. The methodology then proceeds to 
block 124 where the positioning error E r j is determined and recorded. 
Thereafter, the methodology proceeds to decision block 128. 
[0039] In decision block 128 the methodology determines whether the 

value of (i) is equal to zero (0) to thereby determine whether error data has 
been collected for an entire run (r) in the second direction. If the value of (i) 
is not equal to zero (0), the methodology proceeds to block 130 where the 
value of (i) is decremented by one (i.e., i = i - 1). Thereafter, the methodology 
loops back to block 122 so that the controller 16 may increment the tool head 
12 to a location associated with the next axis increment AI M . 
[0040] Returning to decision block 128, if the value of (i) is equal to 

zero (0), the methodology proceeds to decision block 134 where the 
methodology determines whether the value of (r) is equal to a predetermined 
value (R). As the values of E r j are eventually grouped and averaged as will 
be described in detail, below, the value of R is an even number that is greater 
than or equal to four (4). In the example provided, I have assigned R a value 
of six (6) so that three sets of data are taken with the tool head 12 indexed in 
the first direction and three sets of data are taken with the tool head 12 
indexed in the second direction. If the value of (r) is less than (R), the 
methodology proceeds to block 136 where the value of (r) is incremented by 
one (i.e., r = r + 1). Thereafter, the methodology loops back to block 108. 
[0041] Returning to decision block 134, if the value of (r) is equal to (R), 

the methodology proceeds to block 140 where the value of Erev is calculated 
according to the following formula: 
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£ REV — 



\E 2J +E AJ +... + E RJ ) (E U +E 3J +... + E (R _ l)t .) 
(R/2) (R/2) 



0) 



(N + l) 

Alternately, adjusted error values (AE r j) may be substituted for the values of 
E r ,j in the above equation to permit the methodology to utilize compensations 
for cumulative error (CEA), reversal error (REA) and zero positioning error 
(ZPA). The adjusted error value (AE r j) of an associated error value (E r j) 
where the value of (r) is even is calculated as follows: 

AE r j = E r> i - REA + CEA xixL + ZPA (2) 
whereas the adjusted error value (AE r j) of an associated error value (E r j) 
where the value of (r) is odd is calculated as follows: 

AE r j = E r j + CEA x i x L + ZPA. (3) 
When adjusted error values are employed, the resulting value of Erev is useful 
in determining the effect that the various compensations (i.e., cumulative 
error, reversal error and zero positioning error) have on reversal error. 
[0042] The methodology next progresses to block 150 where at least a 
portion of the values of E r j are employed to mathematically develop the line 
EC (shown in Figure 4). Since values of E r j where the value of (r) is even are 
known to include reversal error (which is calculated above), steps must be 
taken to factor this error out so as not to improperly affect the calculation of 
cumulative error. 

[0043] One method by which reversal error may be isolated from the 

calculation of cumulative error is to use only those values of E r j where the 
value of (r) is odd. The slope (m) of the line E c is then provided by the 
equation: 
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2x 
m = 



h (R/2) 



(4) 



(LxN) 

[0044] I prefer, however, to utilize the value of Erev to effectively 
remove the average reversal error from those values of E r ,i where the value of 
(r) is even. In this way, the slope (m) of the line E c is provided by the 
equation: 



2x 
m = 



f jgu + E 3,i + - + E (R-w 1 + Kg^ -Erev) + (E 44 - Erev ) + ... + (E R , - Erev )]} " 



A (5) 



(LxN) 

As in the calculation of reversal error, adjusted error values (AE r ,j) may be 
substituted in the alternative for the values of E r> j in the above equation to 
permit the methodology to utilize compensations for cumulative error (CEA), 
reversal error (REA) and zero positioning error (ZPA). When adjusted error 
values are employed, the resulting value of the slope (m) is useful in 
determining the effect that the various compensations (i.e., cumulative error, 
reversal error and zero positioning error) have on cumulative error. 
[0045] With the value of the slope (m) having been thus calculated, the 
methodology proceeds to block 154 where home error (E H ) is calculated. 
Home error describes the error with which the machine tool 8 positions the 
movable element (e.g., the tool head 12) to a predetermined home position. 
Home error (E H ) is provided by the equation: 

_ {e io -r E 3Q + ... + £^p] + [(£ 2 , 0 -Erev ) + (£ 4>0 - Erev ) + ■■■ + {E R 0 -Erev )]} ^ 

h W) : () 

As noted above for the calculation of reversal error and cumulative error, 
adjusted error values (AE r>i ) may be substituted in the alternative for the 
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values of E r ,j in the above equation to permit the methodology to utilize 
compensations for cumulative error (CEA), reversal error (REA) and zero 
positioning error (ZPA). 

[0046] Thereafter, the methodology proceeds to block 160 where 
revised compensations for cumulative error (RCEA), reversal error (RREA) 
and zero positioning error (RZPE) are calculated as follows: RREA = Erev ; 
RCEA = m; and RZPE = E H . Values for each of the registers 24 and 26 x are 
calculated and preferably downloaded directly from the computer system 38 
to the controller 16. The value of RREA is placed in register 24, while the 
value (v x ) that is placed in a register 26 x is calculated as follows: 

v x = x x d x -RCEA - E H (7) 
[0047] The methodology then proceeds to block 180 where capability 
indices (C pk ) for the machine tool are calculated. As those skilled in the art 
will appreciate, capability indices (C pk ) are employed as a means for 
predicting the ability of a process to consistently produce an output to a 
predefined specification and also function as a problem discovery technique 
for operations running outside specification limits. Generally, the formula for 
Cpk is provided by the following equation: 




(8) 



where Spec is the maximum tolerance, X is the mean value of the data and 
a is the standard deviation of the data. Preferrably, the valued of a is 
calculated through a well known formula for estimating the standard deviation 
of the universe of data points: 

-18- 



HD&P Ref. 7784-000547/COA 
Boeing Ref. No. 00-220B (014666)) 

where S is the sample standard deviation of the data. Although the universe 
from which the data is sampled has an unknown shape, the statistics behind 
the calculation of C pk are sound due to the Central Limit Theorem (which 
states that irrespective of the shape of the distribution of a universe, the 
distribution of averages of subgroups drawn from the universe will tend 
towards normal distribution). 

[0048] In the example provided, the purchasers or manufacturers of a 

machine tool would define a tolerance band for cumulative error, repeatability 
error and position error per unit length and as such, the methodology of the 
present invention calculates three separate C pk indicies: Cumulative C pk , 
Repeat C pk and Position C pk . In the particular example provided, the 
tolerance band for cumulative error is equal to ± 0.001500 inch, the tolerance 
band for repeatability error is equal to ± 0.000500 inch and the tolerance band 
for positioning error per unit length is equal to ± 0.000700 inch / 12 inches, 
which are relatively tight or close tolerances as the machine tool 8 is 
employed in the fabrication of critical aerospace tooling, components and 
such. Those skilled in the art will appreciate that the tolerance bands used for 
another CNC machine or CMM may vary from those that are provided here 
depending on the type of machine tool and the nature of the work that the 
machine tool is to perform. 

[0049] As the registers 24 and 26 x have been programmed in a manner 
that compensates for assignable causes of error to thereby center the 
process, the values for (Spec) that are used in calculating C pk indices are 
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simply one-half of the associated tolerance band. Accordingly, the Spec 
values for cumulative (Spec C uM), repeatability (Spec R E P ) and position error per 
unit length (Spec PO s) are 0.001500 inch, 0.000500 inch and 0.000700 inch / 
12 inches, respectively. 

[0050] Relying on the above formula for Cumulative C pk is calculated as 

follows: 



Cumulative C pk = 



S P ec cuM -\Xcm 



3x<r 



(10) 



CUM 



where: 



CUM — ■ 



N 


R 


z 




1=0 


R 



(N + l) 



(11) 



Scum ~ 



i-0 



R £ 

/ J — — X CUM 



~\2 



r=l 



N 



(13) 



and 



a cuM ~ $cum X a/ 1 + 



1 



(14) 



It should be noted that the total number of data samples is (N+1) so that the 
formula for Scum and a CUM nonetheless properly calculates standard deviation 
&est for cumulative positioning. As noted above for the calculation of reversal 
error, cumulative error and home error, adjusted error values (AE r ,j) may be 
substituted in the alternative for the values of E r>j in the two immediately 
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preceding equations to permit the methodology to utilize compensations for 
cumulative error (CEA), reversal error (REA) and zero positioning error (ZPA). 
[0051] Similarly, Repeat C pk may be calculated as follows: 



where: 



and 



Repeat C pk = 



S P eC REP - \XrEP\ 



3xa 



REP 



(15) 



i=0 



X REP — ' 



r=2 



(/?-!) 



(16) 




(17) 



a REP ~ SrEP X i| ^ + 



1 



N+l 



(18) 



As noted above, adjusted error values (AE r ,j) may be substituted in the 
alternative for the values of E r j in the two immediately preceding equations to 
permit the methodology to utilize compensations for cumulative error (CEA), 
reversal error (REA) and zero positioning error (ZPA). 
[0052] Also similarly, Position C pk may be calculated as follows: 

Spec^ - Xpos 

Position C pk = ! 1 (19) 

3 x <t pos 

where: 
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POS 



= 1 

1=0 



R F -h F 

r=l 2 



(20) 




(21) 



and 



°>OS — ^ POS 



xJfi + -L- 



(22) 



As noted above, adjusted error values (AE rJ ) may be substituted in the 
alternative for the values of E r j in the two immediately preceding equations to 
permit the methodology to utilize compensations for cumulative error (CEA), 
reversal error (REA) and zero positioning error (ZPA). 

[0053] The methodology preferably also calculates "safe" tolerances for 
cumulative error, repeatability and position per unit length. A workpiece 
having tolerances that exceed their corresponding safe tolerance may be 
accurately and repeatably fabricated by the machine tool 8. Generally 
speaking, a safe tolerance is calculated as follows: 

ST = ±(3xcj + \x\) (23) 
where associated values of crand Xare employed for a given error {e.g., 



'CUM 



and Xcum are employed to calculate a safe tolerance for cumulative 



error (ST CU m)}. Thereafter, the methodology progresses to bubble 200 and 
terminates. 
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[0054] Due to the statistics behind their calculation, the values of 

(Erev), (m) and (E H ), Capability Indices (Cpks) and Safe Tolerances (ST) are 
relatively robust in that they may be readily verified in a subsequent capability 
analysis (i.e., a second capability analysis conducted immediately after a first 
capability analysis). Accordingly, the values of (E RE v), (m) and (Eh), 
Capability Indices (Cpks) and Safe Tolerances (ST) are important from a 
historical perspective, as they document the performance of an associated 
axis of the machine tool. In this regard, performance of a capability analysis 
at regular intervals using the methodology of the present invention will inform 
the mechanic of changes in the machine tool that affect its performance but 
which are otherwise extremely difficult (or even impossible) to identify. Stated 
another way, by conducting a regularly scheduled capability analysis using 
the methodology of the present invention and comparing the values of (£W), 
(m) and (Eh), Capability Indices (Cpks) and Safe Tolerances (ST) to those 
obtained in one or more previously conducted capability analysis (using the 
methodology of the present invention), a mechanic is able to identify very 
subtle changes in the axis or axes of the machine tool which affect 
performance but which are otherwise difficult or impossible to trace back to a 
particular axis or axes. 

[0055] While the methodology of the present invention has been 
described thus far as calculating the values of (Erev), (m) and (E H ) directly, 
those skilled in the art will appreciate that the invention, in its broader aspects, 
may be practiced somewhat differently. For example, the methodology may 
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alternately calculate the values of (Erev) and (m) and EH in an iterative 
manner as is illustrated in Figure 6. 

[0056] In this embodiment, the methodology is generally similar to the 

methodology of Figure 5 except for the addition of decision block 170 and 
block 174 and the use of adjusted error values (AE r ,j) is mandatory rather than 
optional. When the methodology completes block 160, it progresses to 
decision block 170 to determine whether the values of Erev , m and E H have 
been optimized or reduced to a point where they are less than or equal to a 
predetermined threshold. In my experience, I have found that the values of 
Erev , m and E H can be reduced to about zero (0), but the predetermined 
threshold may be higher if desired. Furthermore, the predetermined 

thresholds associated with the optimization of the values of ~Erev , m and E H 
may be dissimilar. 

[0057] Accordingly, if the values of E RE v , m and E H are not less than or 
equal to their respective predetermined thresholds, the methodology 
progresses to block 174 where the values of CEA, REA and ZPA are 
recalculated as follows: 

CEA = CEA + RCEA (24) 
REA = REA + RREA (25) 
ZPA = ZPA + RZPA. (26) 
The methodology then loops back to block 1 40 so that Erev , m and E H may 
be recalculated with the recalculated values of CEA, REA and ZPA. 
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[0058] Returning to decision block 170, if the values of Erev , m and E H 

are less than or equal to their respective predetermined thresholds, the 
methodology progresses to block 180 and proceeds as described above. 
[0059] While the invention has been described in the specification and 
illustrated in the drawings with reference to a preferred embodiment, it will be 
understood by those skilled in the art that various changes may be made and 
equivalents may be substituted for elements thereof without departing from the 
scope of the invention as defined in the claims. In addition, many modifications 
may be made to adapt a particular situation or material to the teachings of the 
invention without departing from the essential scope thereof. Therefore, it is 
intended that the invention not be limited to the particular embodiment illustrated 
by the drawings and described in the specification as the best mode presently 
contemplated for carrying out this invention, but that the invention will include 
any embodiments falling within the foregoing description and the appended 
claims. 
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